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Trapping ultracold atoms in optical lattices enabled numerous breakthroughs uniting several dis¬ 
ciplines. Coupling these systems to quantized light leads to a plethora of new phenomena and has 
opened up a new field of study. Here we introduce a physically novel source of competition in a 
many-body strongly correlated system: We prove that quantum backaction of global measurement 
is able to efficiently compete with intrinsic short-range dynamics of an atomic system. The compe¬ 
tition becomes possible due to the ability to change the spatial profile of a global measurement at 
a microscopic scale comparable to the lattice period without the need of single site addressing. In 
coherence with a general physical concept, where new competitions typically lead to new phenom¬ 
ena, we demonstrate novel nontrivial dynamical effects such as large-scale multimode oscillations, 
long-range entanglement and correlated tunneling, as well as selective suppression and enhancement 
of dynamical processes beyond the projective limit of the quantum Zeno effect. We demonstrate 
both the break-up and protection of strongly interacting fermion pairs by measurement. Such a 
quantum optical approach introduces into many-body physics novel processes, objects, and meth¬ 
ods of quantum engineering, including the design of many-body entangled environments for open 
systems. 


I. INTRODUCTION 

Ultracold gases trapped in optical lattices created by 
classical laser beams represent a successful interdisci¬ 
plinary field: Atomic systems allow quantum simulations 
of phenomena predicted in condensed matter and parti¬ 
cle physics, and find applications in quantum informa¬ 
tion processing [T]. Coupling quantum gases to quan¬ 
tized light (cf. Refs. [1113] and references therein) broad¬ 
ens the field even further with opportunities to obtain 





FIG. 1. Setup. Atoms in an optical lattice are probed by 
a coherent light beam, and the light scattered at a particular 
angle is enhanced and collected by a leaky cavity. The pho¬ 
tons escaping the cavity are detected, perturbing the atomic 
evolution via measurement backaction. This process com¬ 
petes with the usual (Bose-)Hubbard dynamics and leads to 
the novel phenomena described in the text. 
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novel quantum phases and light-matter entanglement 01 - 
|7]. Here we propose how the phenomenon of quantum 
measurement backaction can be made efficiently com¬ 
pete with standard short-range processes in optical lat¬ 
tices (tunneling and on-site interaction), with resulting 
novel nontrivial effects, which do not require additional 
quantum control |S]. This fits well into a general phys¬ 
ical concept, where a novel competition in many-body 
systems typically leads to new phenomenon. The col¬ 
lapse of the atomic wave function by detecting scattered 
light reflects one of the most fundamental manifestations 
of quantum mechanics - measurement backaction - 
here due to light-matter entanglement. We go beyond 
recent quantum nondemolition (QND) approaches |10I - 
I18| . where either many-body dynamics or measurement 
backaction did not play any role. We demonstrate that 
the key mechanism enabling us to construct efficient com¬ 
petition between the global backaction and short-range 
processes is the possibility to spatially structure the mea¬ 
surement at a microscopic scale comparable to the lat¬ 
tice period without the need for single site resolution. 
By carefully choosing the geometry of light we can select 
which subspace of atomic dynamics is affected by mea¬ 
surement. We consider a flexible setup where the global 
light scattering can be engineered, allowing us to sup¬ 
press or enhance specific dynamical processes, realizing 
a spatially nonlocal quantum Zeno effect. The measure¬ 
ment process generates spatial modes m of matter fields 
that can be considered as designed systems and reser¬ 
voirs, opening the possibility of controlling dissipations 
in ultracold atomic systems, without resorting to atom 
losses and collisions which are difficult to manipulate. 
The continuous measurement of the light field introduces 
a controllable decoherence channel into the many-body 
dynamics. Global light scattering from multiple lattice 
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sites creates nontrival, spatially nonlocal coupling to the 
environment, which is impossible to obtain with local 
interactions |20H23| . Such a quantum optical approach 
can broaden the field even further, allowing quantum 
simulation of models unobtainable using classical light, 
and the design of novel systems beyond condensed mat¬ 
ter analogs. For example, both designed systems and 
reservoirs are represented by many-body strongly cor¬ 
related systems with internal long-range entanglement. 
This will raise novel type of questions and stimulate fur¬ 
ther research on dynamics and states in physics of open 
systems and quantum engineering. Additionally, many- 
body processes in dissipative systems have recently at¬ 
tracted attention in biological systems |24| . 

In contrast to other approaches of studying ground or 
steady state properties, here we focus on the dynamics of 
single quantum trajectories and find effects that are not 
visible in the ground and steady states. For bosons, the 
measurement induces dynamical macroscopic superposi¬ 
tions (multimode Schrodinger cat states). The resulting 
states have both oscillating density-density correlations 
with nontrivial spatial periods and long-range coherence, 
thus having properties of the supersolid state [25], but 
in the essentially dynamical version. For strongly in¬ 
teracting fermions, we demonstrate how to destroy and 
protect fermion pairs, clearly showing the interplay be¬ 
tween many-body dynamics of fermions and the quantum 
measurement backaction. Even in the strong measure¬ 
ment regime, where the backaction dominates the sys¬ 
tem evolution and projects it onto a single eigenstate, 
we find that long-range correlated tunneling emerges, 
leading to dynamical generation of long-range entangle¬ 
ment between distant sites (which can exist even in one¬ 
dimensional systems). The setup we consider is analo¬ 
gous to the famous cavity QED experiments |3|, where 
the quantum light states in a cavity were probed and pro¬ 
jected by measuring atoms. Here, light and matter are 
reversed and we probe matter-fields with light. A strik¬ 
ing advancement is that the number of quantum matter 
waves (sites with trapped atoms) can be easily scaled 
from few to thousands by changing the number of illumi¬ 
nated sites, while scaling cavities is a challenge [^. This 
work, together with recent experiments where our predic¬ 
tions can be tested (Bose-Einstein Condensates trapped 
inside a cavity |25fl27| . and recently a lattice in a cavity 
|28l [2^ 1. will help to develop a deeper understanding of 
the effect of quantum measurement on the dynamics of 
strongly correlated many-body systems and uncover the 
constructive and active role of global backaction in sys¬ 
tems with short-range interactions, thus closing the gap 
between quantum optics and quantum gases. 


II. MODEL FOR SPATIALLY STRUCTURED 
GLOBAL MEASUREMENT 

We consider off-resonant light scattering from N atoms 
trapped in an optical lattice with period d and L sites |5] 


(see Appendix A and Fig. [^. The light scattered at a 
particular angle can be selected and enhanced by a cav¬ 
ity |30fl5^ with decay rate k. Similar to classical optics, 
the light amplitude is given by a sum of scatterings from 
all atoms with coefficients dependent on their positions: 
a = C{D + B), where a is the photon annihilation opera¬ 
tor, C is the Rayleigh scattering coefficient (see Appendix 
A) and 

L L 

D = Jjjfi'j, B = Jijh\hj, (1) 

i=l (ij> 

where bj and fij = h^-bj are the atomic annihilation and 
number operators at site j sums over neighboring 

sites) and Jij are given by 

Ji] = jw{v- r,)u*^^{r)uin{r)w{r - rj) dr, (2) 

where w{r) are the localized Wannier functions and 
Win, out (r) are the mode functions of incoming and out¬ 
going light respectively (e.g., ui{r) = exp(jk;-r) for trav¬ 
eling and ui(r) = cos(ki-r) for standing waves with wave 
vectors k;). 

In Eq. ([^ D describes scattering from the on¬ 
site densities, while B that from the inter-site coher¬ 
ence terms [33]. For well-localized atoms, the second 
term is usually neglected, and a = CD with Jjj = 
u*ut(rj)win(rj). For spin-^ fermions we use two light 
polarizations ax,y that couple differently to two spin 
densities n^j, allowing measurement of their lin¬ 
ear combinations, e.g., = CD^ = 

Gy = CDy = Jjjihj, where pj = fi^j + n^j and 

ifij = rt/fj — fiij are the mean density and magnetisa¬ 
tion. This property has recently been used to investigate 
spin-spin correlations in Fermi gases |34l 135] . 

We focus on a single run of a continuous measure¬ 
ment experiment using the quantum trajectories tech¬ 
nique [33] (see Appendix B). The evolution is determined 
by a stochastic process described by quantum jumps (the 
jump operator c = y/^a is applied to the state when 
a photodetection occurs) and non-Hermitian evolution 
with the Hamiltonian 

= Ho - ihJ c/2 (3) 

between jumps, where Hq is the usual (Bose-)Hubbard 
Hamiltonian. Importantly, the measurement introduces 
a new energy and time scale 7 = IC^k, which competes 
with the two other standard scales responsible for unitary 
dynamics of closed systems (tunneling J and on-site in¬ 
teraction U). If atoms scatter light independently, inde¬ 
pendent jump operators Cj would be applied to each site, 
projecting the atomic system to a state, where the long- 
range coherence degrades Ezj. This is a typical scenario 
for spontaneous emission I3ZIISH], or rather analogous 
local and fixed-range |43l IHj addressing 
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and interactions. Additionally, if the light is scattered 
without a cavity, i.e. in uncontrolled directions, we lose 
the ability to choose the measurement operator and the 
jump operator applied would depend on the direction of 
the detected photon. In contrast, here we consider global 
coherent scattering, where the single global jump oper¬ 
ator c is given by the sum over all sites, and the local 
coefficients Jjj Q responsible for the atom-environment 
coupling (via the light mode a) can be engineered by 
optical geometry. Thus, atoms are coupled to the en¬ 
vironment globally, and atoms that scatter light with 
the same phase are indistinguishable to light scattering 
(i.e. there is no “which-path information”). As a strik¬ 
ing consequence, the long-range quantum superpositions 
are strongly preserved in the final projected states, and 
the system splits into several spatial modes m, where 
all atoms belonging to the same mode are indistinguish¬ 
able, while being distinguishable from atoms belonging 
to different modes. 

We engineer the atom-environment coupling coeffi¬ 
cients Jjj using standing or traveling waves at different 
angles to the lattice. A key mechanism, which will allow 
us to construct the effective competition between global 
measurement and local processes, is the ability to modify 
these couplings at very short microscopic distances. This 
is in striking contrast to typical scenarios of Dicke and 
Lipkin-Meshkov-Glick models |3^, where the coupling is 
global, but rather homogeneous in space. If both probe 
and scattered light are standing waves crossed at such an¬ 
gles to the lattice that projection ki„-r is equal to kout’T 
and shifted such that all even sites are positioned at the 
nodes (do not scatter light), one gets Jjj = 1 for odd 
and Jjj = 0 for even sites. Thus we measure the number 
of atoms at odd sites only (the jump operator is propor¬ 
tional to IVodd)) introducing two modes, which scatter 
light differently: odd and even sites. The coefficients 
Jjj = (—1)-^ are designed by crossing light waves at 90° 
such that atoms at neighboring sites scatter light with tt 
phase difference |TU1 ITFl . giving D = TVeven — .^odd, 

introducing the same modes, but with different coherence 
between them. Moreover, using travelling waves crossed 
at the angle such that each i?-th site is indistinguish¬ 
able ((kin — kout)'rj = 27rj /R), introduces R modes with 
macroscopic atom numbers Nf. D = |19| . 

Here two (odd- and even-site modes) appear for R = 2. 
Therefore, we reduce the jump (measurement) operator 
from being a sum of numerous microscopic contributions 
from individual sites to the sum of smaller number of 
macroscopically occupied modes with a very nontrivial 
spatial overlap between them. 

Our results are applicable in any number of dimen¬ 
sions and for an arbitrary system size, but for simplic¬ 
ity and clarity our results are often presented in one di¬ 
mension. The dimensionality of the system affects key 
properties of the ground state such as quantum criti¬ 
cal points and decay of atomic correlations. However, 
in this paper the ground state represents only a realistic 
initial condition and the long-range correlations are dom- 



FIG. 2. Large oscillations between the measurement-induced 
spatial modes resulting from the competition between tnnnel- 
ing and weak-measurement backaction. The plots show sin¬ 
gle quantum trajectories, (a)-(d), Atom number distributions 
p{Ni) in one of the modes, which show various number of well- 
squeezed components, reflecting the creation of macroscopic 
superposition states depending on the measurement configu¬ 
ration {U/ J = 0, 'y/ J — 0.01, L = N, initial states: superfluid 
for bosons, Fermi sea for fermions), (a), Measurement of the 
atom number at odd sites NoAd creates one strongly oscillat¬ 
ing component in p(Aodd) {N = 100 bosons, Jjj = 1 if j is 
odd and 0 otherwise), (b). Measurement of (iVodd — fVeven)^ 
introduces R = 2 modes and preserves the superposition of 
positive and negative atom number differences in p{Nodd) 
{N = 100 bosons, Jjj = (—l)-’^^). (c). Measurement for 

R = 3 modes (see text) preserves three components in p{Ni) 
{N = 108 bosons, Jjj = (d). Measurement of fVodd 

for fermions leads to oscillations in p(A^odd), though not as 
clearly defined as for bosons because of Pauli blocking (L = 8 
sites, Nj = Nj = 4 fermions, Jjj = 1 if j is odd and 0 other¬ 
wise). Simulations for ID lattice. 


inated by the effect of the nonlocal nature of the mea¬ 
surement. Therefore, the results that we will present do 
not strictly depend on the number of spatial dimensions 
but on the geometry and symmetry of the measurement 
scheme. The only exception is represented by the behav¬ 
ior of ID fermions: This differs from any other dimension 
since atoms with the same spin can not pass each other. 

In the following, we will show how globally designed 
measurement backaction introduces spatially long-range 
interactions of the modes, and demonstrate novel effects 
resulting from the competition of mode dynamics with 
standard local processes in a many-body system. 
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III. LARGE-SCALE DYNAMICS DUE TO 
WEAK MEASUREMENT 

We start with non-interacting bosons {U/J = 0), and 
demonstrate that the competition between tunneling and 
global measurement strongly affects the dynamics of the 
atomic system. The weak measurement (7 <C J) is un¬ 
able to freeze the atom numbers projecting the atomic 
state and quantum Zeno dynamics |491155] cannot be es¬ 
tablished. In contrast, the measurement leads to giant 
oscillations of particle number between the modes. Fig¬ 
ures [ 2 ^-c illustrate the atom number distributions in one 
of the modes for R = 2 (Nodd) and R = 3 (Ni). Without 
continuous monitoring, these distributions would spread 
out significantly and oscillate with an amplitude propor¬ 
tional to the the initial imbalance, i.e. tiny oscillations for 
a tiny initial imbalance. In contrast, here we observe (i) 
full exchange of atoms between the modes independent 
of the initial state and even in absence of initial imbal¬ 
ance, (ii) the distributions consist of a small number of 
well-defined components, and (iii) these components are 
squeezed even by weak measurement. Depending on the 
quantities addressed by the measurement, the state of 
the system has a multi-component structure which is a 
consequence of the photon number (intensity) a^a not 
being sensitive to the light phase. In other words, the 
measurement does not distinguish between all permu¬ 
tations of mode occupations that scatter light with the 
same intensity. The number of components of the atomic 
state, i.e. the degeneracy of a^a, can be computed from 
the eigenvalues of noting that they 

can be represented as the sum of vectors on the complex 
plane with phases that are integer multiples of 27r/i?: 
]\rig* 2 -n-/fl, Since the sum of these vec¬ 

tors is invariant under rotations by 2^1/R, I S Z and 
reflection in the real axis, the state of the system is 2 - 
fold degenerate for i? = 2 and 2 ii-fold degenerate for 
R > 2. Figure shows the superposition of two states 
with positive and negative Aodd ~ -^evem while Fig. H 
illustrates the superposition in the three-mode case. 

One can get a physical insight into the origin of oscil¬ 
lations by constructing the following model. For macro- 
scopically occupied modes of non-interacting atoms ini¬ 
tially in the superfiuid state, the i?-modes problem can 
be treated analytically reducing it to R effective sites. In 
particular, for i? = 2 (effective double-well |54l [hH] !. we 
can write the atomic state as 

N 

\i:)=Y,qi\l,N-l). (4) 

1=0 

where the ket \l,N — 1) represents a superfiuid with I 
atoms in the first and N — I atoms in the second spatial 
mode. In the limit A 1 we can describe the evolution 
of the system using continuous variables |55| and define 
the wave function = l/N) = y/Nqi. Introducing 
the relative population imbalance between the two wells 
z = 2 a: — 1 and starting from the superfiuid state, the 


solution of the Schrodinger equation with Hamiltonian 
iJeff is 


' 0 (z, t) oc exp 


ia{t) -\- 


izc{t) 

W(i) 


262 (t) 


{Z-Zd{t)f 

262 (t) J ’ 

(5) 


where h(t) is the width of atomic distribution, Zdit) is 
the average population imbalance and c{t) and are 
real phases while a{t) is complex and takes into ac¬ 
count the decay of norm of the wave function due to 
the non-Hermitian term present in Hes- Expanding 
the Hamiltonian in powers of 1/A up to second order 
we obtain a system of coupled differential equations for 
a(t), b(t), c{t), zo(t), and which, in the weak mea¬ 
surement limit, can be linearized, leading to 

X IV *7 X 

zo{t) = -e“T“‘ [c(0) sin(2Jt) -|- 2zo(0) cos(2Jt)]. ( 6 ) 

This expression describes the evolution of the popula¬ 
tion imbalance between two quantum jumps: The atoms 
oscillate between the two spatial modes with decreasing 
amplitude and tend to restore a balanced distribution. 
The quantum jumps strongly affect the dynamics of the 
system via measurement backaction, driving the oscilla¬ 
tions to their possible maximum amplitude, i.e., all the 
atoms oscillate between the two spatial modes. We ex¬ 
plain this effect noting that (i) the average time between 
two jumps (~ 2 /A^ 7 ) is much smaller than the damping 
time in (|^, (ii) the probability of a jump depends on the 
imbalance itself pjump cx; (1 -f zq)^ and (iii) each jump 
tends to increase the population imbalance. In order to 
calculate how the photodetection events affect the evo¬ 
lution of zq) we compute the effect of average jump rate 
as 6zo/St = Pjump ^zo/St, where Azq is the change in zq 
due to the jump itself. Solving this equation leads to 

2 ^ 0 ,jumps(l) = + (1 + ZQ{0))e^^^. (7) 

Therefore, the amplitude of oscillations of zq increase 
since the difference between the exponents in (|^ and 
([^ is positive and the system leaves the stable point 
Zq = z'o = 0. Note that the dynamics described here 
is a feature of single trajectories only. This is in con¬ 
trast to averaging over many runs, which corresponds to 
the master equation solution, which masks these effects 
completely. This happens because the oscillation phase 
changes from realization to realization as is known from 
works involving single and multiple measurements m- 

The competition between measurement and atomic dy¬ 
namics allows realization of multicomponent macroscopic 
superpositions (Schrodinger cat or NOON states), which 
are useful in quantum metrology and information. Such 
multimode superpositions are a purely quantum effect. 
The multimode dynamics may recover the semiclassical 
character m, when the number of modes is reduced to 
1. The method we propose does not require external 
control jSl UHl |59] for preparing these states: By contin¬ 
uously monitoring the light intensity it is possible to de¬ 
termine when the splitting in the components reaches its 
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maximum value (corresponding to the maximal macro- 
scopicity m) and further oscillatory dynamics can be 
stopped by ramping up the lattice depth. Note that for 
R > 1 spatial modes each photocount changes the phase 
difference between the various components of the atomic 
state, making it fragile to photon losses. However, the 
measurement setup can be modified to make these states 
more robust m- In addition, the example in Fig. ga) 
consists of only one component (here A^odd is measured 
directly) and is therefore insensitive to decoherence due 
to photon losses. 

It is interesting to note that the oscillating state we 
have described shows spatial periodicity in the density- 
density correlation function depending on the measure¬ 
ment configuration. For R modes the spatial period is 
Rd {d is the lattice period). Moreover, the state also has 
long-range coherence between distant sites. Therefore, it 
has properties of the supersolid state |25] but with a non¬ 
trivial period, and it exists in the essentially dynamical 
version. 

In contrast to bosons, dynamics for two modes of non¬ 
interacting fermions does not show well-defined oscilla¬ 
tions (Fig. il) due to Pauli exclusion. However, while 
the initial ground state is a product of f and wave func¬ 
tions (Slater determinants), the measurement introduces 
an effective interaction between two spin components and 
the state becomes entangled by measurement [51]. The 
dynamics of a one-dimensional fermionic system strongly 
depends on the spatial profile of the measurement oper¬ 
ator. For example, if only the central part of the lattice 
is illuminated (diffraction maximum) the atom transfer 
between the modes induced by the measurement is sup¬ 
pressed as the presence of an atom at the edges of the 
illuminated area completely forbids the atomic tunneling 
between the two modes, greatly decreasing the fluctua¬ 
tions of the measurement operator. 

Carefully choosing geometry, one can suppress the on¬ 
site contribution to light scattering and effectively con¬ 
centrate light between the sites, thus in-situ measuring 
the matter-field interference [33]. In this case, a = 

CB, and the coefficients Jij for i ^ j from Eq. Q can be 
engineered |33| . For = 1, the jump operator is propor¬ 
tional to the kinetic energy Ek = —2hJ J2k cos{ka) 

and tends to freeze the system in eigenstates of the non¬ 
interacting Hamiltonian. The measurement projects to a 
superposition of two states with different kinetic energies: 
a superposition of matter waves propagating with differ¬ 
ent momenta. The measurement does not distinguish 
between these two states, because in a lattice, the two 
momentum states |fc) and | 7 r/a — fc) interfere in the same 
way, but with opposite phase in between the lattice sites. 
The measurement freezes dynamics for any y/J, since 
the jump operator and Ffeff have the same eigenstates. 
As a result of the detection, the atoms quickly spread 
across the lattice, and the density distribution becomes 
uncertain (Fig. clearly illustrating the quantum un¬ 
certainty relation between the number- and phase-related 


variables {fii and b\bi+i). Note that, in absence of mea¬ 
surement, such a distribution presents a periodic spread 
and revival due to coherent tunneling and the system 
does not reach a steady state (Fig. Therefore, engi¬ 
neering Jij can lead to the measurement-based prepara¬ 
tion of peculiar multicomponent momentum (or Bloch) 
states. 


IV. COMPETITION BETWEEN 
INTERACTION, TUNNELING AND 
MEASUREMENT 

As we turn on the inter-atomic interactions, U/J 7 ^ 0, 
the atomic dynamics changes as the measurement com¬ 
petes with both the tunneling and this on-site interaction. 
One approach is to study the ground or steady state of 
the system in order to map a quantum or dissipative 
phase diagram. This is beyond the scope of this arti¬ 
cle, because here, we adopt a quantum optical approach 
with the focus on the conditional dynamics of a quantum 
trajectory corresponding to a single experimental realiza¬ 
tion. The resulting evolution does not necessarily reach 
a steady state and can occur far from the ground state 
of the system. Again, each quantum trajectory evolves 
differently as the detection process is determined stochas¬ 
tically and even states with similar expectation values of 
D can have minimal overlap. However, even though each 
trajectory is different they all have one feature in com¬ 
mon: The uncertainty in the measured operator, D, is 
only a function of the Hamiltonian parameters, 7 , J, and 
U. Therefore, we average its variance over many real¬ 
izations ((CT|))traj) as this quantity effectively describes 
the squeezing of the atomic distribution due to measure¬ 
ment. Importantly, it is not possible to access this quan¬ 
tity using the master equation solution: The uncertainty 
in the final state is very large and it completely hides 
any information on the spread of a single trajectory. In 
other words, the master equation addresses the variance 
of the average value of D over the trajectories ensemble 
((-D^)traj — (-O)traj) Squeezing of a single tra¬ 

jectory conditioned to the measurement outcome. This 
again highlights the fact that interesting physics happens 
only at the single trajectory level. In this section we show 
results for the measurement oi D = which is ro¬ 

bust to photon losses. Specifically, we compute the width 
of the atomic distribution (((T|))traj) in the limit t —>■ oo, 
when its value does not change significantly in time even 
if the atomic imbalance is not constant. Moreover, we use 
the ground state of the system as initial state since this 
is a realistic starting point and a reference for explaining 
the measurement-induced dynamics. 

For bosons (Fig. the number fluctuations cr|j cal¬ 
culated in the ground state decrease monotonically for 
increasing U, reflecting the superfluid - Mott insulator 
quantum phase transition. The measured state on the 
other hand behaves very differently and (cr|j)traj varies 
non-monotonically. For weak interaction, the fluctua- 
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FIG. 3. Evolution of the on-site atomic density, while measur¬ 
ing the matter-field coherence between the sites B. The plots 
show single quantum trajectories, (a), Atoms, all initially at 
the edge site, are quickly spread across the whole lattice lead¬ 
ing to the large uncertainty in the atom number, while the 
matter-phase related variable is defined (projected) by the 
measurement (bosons, N = L = 6, U/J = 0, q/J = 0.1, 
Jjj ~ !)• (b), Atomic density spread and revival due to 

coherent tunnelling in the absence of measurement (bosons, 
N = L = 6, U/J = 0, 7 /J = 0, Jjj = 1). Simulations for ID 
lattice. 


tions are strongly squeezed below those of the ground 
state; then they quickly increase, reach their maximum 
and subsequently decrease as the interaction becomes 
stronger. We explain this effect looking at the dynam¬ 
ics of single trajectories (cf. Fig. [^. For small values 
of U/J, the population imbalance between odd and even 
sites oscillates and its uncertainty is squeezed by the mea¬ 
surement as described in Section III. However, when such 
oscillations reach maximal amplitude the local atomic re¬ 
pulsion spreads the atomic distribution and prevents the 
formation of states with large atom number in one of the 
two modes (Fig. Since the interaction is a nonlinear 
term in the atomic dynamics, states with different im¬ 
balance oscillate with different frequencies and the mea¬ 
surement is not able to squeeze the fluctuations of iVodd 
as efficiently as in the non-interacting case. This behav¬ 
ior in the weak measurement and weak interaction limit 
is in contrast with what the effective two-sites model in¬ 
troduced in the previous section predicts (in a double 
well fluctuations are monotonically decreasing with U). 
However, the solution to this model is only valid in the 
^ 1 limit which suggests that this difference is due to 
the fact that in a double well with a large occupancy the 
population transfer is sequential, i.e. atoms transfer one 
by one from one well to the other, increasing the total 
interaction energy, which goes up as (n^), very steadily. 
On the other hand, in a lattice a collective excitation of 
a single atom from each site in one mode to another site 
in the other mode will increase the energy by KU, where 
K is the number sites in one mode. This is an increase 
by a factor of K compared to a single particle-hole exci¬ 
tation pair. Therefore, a lattice and a global, long-range 
measurement scheme are necessary to observe such a col¬ 
lective transfer of atoms. 


FIG. 4. Atom number fluctuations demonstrating the com¬ 
petition of global measurement with local interaction and tun¬ 
neling. Number variances are averaged over many trajecto¬ 
ries. Error bars are too small to be shown (~ 1%), which 
emphasizes the fundamental nature of the squeezing, (a), 
Bose-Hubbard model with repulsive interaction. The fluctua¬ 
tions of the atom number at odd sites Aodd in the ground 
state without a measurement (thin solid line) decrease as 
U/J increases, reflecting the transition between the super¬ 
fluid and Mott insulator phases. For the weak measurement, 
(o'£))traj is squeezed below the ground state value, but then 
increases and reaches its maximum as the atom repulsion pre¬ 
vents oscillations and makes the squeezing less effective. In 
the strong interacting limit, the Mott insulator state is de¬ 
stroyed and the fluctuations are larger than in the ground 
state. (100 trajectories, N = L — G, Jjj = 1 if i is odd 
and 0 otherwise.) (b),(c). Fermionic Hubbard model with 
attractive interaction; fluctuations of the total atom number 
at odd sites = iV^odd + iV 4 .odd (b) and of the magnetiza¬ 
tion at odd sites Dy = Modd = iV^odd — iV^odd (c). Without 
measurement, the interaction favors formation of doubly oc¬ 
cupied sites so the density fluctuations in the ground state are 
increasing while the magnetization ones are decreasing. The 
measurement creates singly occupied sites decreasing the den¬ 
sity fluctuations and increasing the magnetization ones, which 
manifests the break-up of fermion pairs by measurement. The 
measurement-based protection of fermion pairs is shown in the 
next figure. (100 trajectories, L = 8, Nj — = 4, Jjj = 1 if 

j is odd and 0 otherwise.) Simulations for ID lattice. 


For weak measurement, but in the strongly interacting 
limit, we note that the measurement leads to a significant 
increase in fluctuations compared to the ground state. 
Both, the local interaction and measurement squeeze 
fluctuations, but as the measurement destroys the Mott 
insulator, the fluctuations are larger than in the ground 
state. Using first-order perturbation theory the ground 
state of the system is 


\^j/u) 


1 + 77 


u 


(iJ) 




( 8 ) 


where j'l'o) is the Mott insulator state and the second 
term represents a uniform distribution of particle-hole 
excitation pairs across the lattice. The action of a single 
photocount will amplify the present excitations increas- 
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ing the fluctuations in the system. In fact, consecutive 
detections lead to an exponential growth of these excita¬ 
tions as for if ^ 1 and unit Ailing, and the atomic state 
after m quantum jumps becomes j/u) oc |'I'j/; 7 ) + 
|$m) where 

i odd 

(9) 

In the weak measurement regime the effect of the non- 
Hermitian decay is negligible compared to the local 
atomic dynamics combined with the quantum jumps and 
so there is minimal dissipation occurring. Therefore, be¬ 
cause of the exponential growth of the excitations, even 
a small number of photons arriving in succession can de¬ 
stroy the Mott insulator state very quickly. This will 
always happen given sufficient time, and provided there 
are finite fluctuations present in the initial state and so 
this will happen at any value of U/J, except for J = 0 
when the ground state becomes a single Fock state with 
no particle-hole excitations. 

In the strong measurement regime (7 ^ J) the mea¬ 
surement becomes more significant than the local dynam¬ 
ics and the system will freeze the state in the measure¬ 
ment operator eigenstates. In this case, the squeezing will 
always be better than in the ground state, because mea¬ 
surement and on-site interaction cooperate in suppressing 
fluctuations. For low interaction strengths this should 
be obvious from the fact that in a superfluid ground 
state the atoms are spread out over the entire lattice and 
thus the uncertainty in atom number is large whereas 
measurement eigenstates have a well-defined occupation 
number. However, the strongly interacting regime is 
much less evident, especially since we have demonstrated 
how sensitive the Mott insulating state is to the quantum 
jumps when the measurement is weak. 

To understand the strongly interacting case we will 
again use first-order perturbation theory and consider a 
postselected {&£>) = 0 trajectory. This corresponds 
to a state that scatters no photons and so the non- 
Hermitian correction to the Hamiltonian is sufficient to 
understand the measurement. Since squeezing depends 
on the measurement strength and is common to all pos¬ 
sible trajectories we can gain some insight by consider¬ 
ing this specific case. However, we will now consider 
D = AN = Nodd — Neven as this measurement also has 
only two modes, R = 2, but its D) = 0 trajectory cor¬ 
responds to the Mott insulating ground state. According 
to perturbation theory the modified ground state is now 



FIG. 5. Conditional dynamics of the atom-number distribu¬ 
tions at odd sites illustrating competition of the global mea¬ 
surement with local interaction and tunneling (single quan¬ 
tum trajectories, initial states are the ground states), (a), 
Weakly interacting bosons: the on-site repulsion prevents the 
formation of well-defined oscillations in the population of the 
mode. As states with different imbalance evolve with dif¬ 
ferent frequencies, the squeezing due to the measurement is 
not as efficient as one observed in the non-interacting case 
{N = L = 6, U/ J = 1, 7 /J = 0.1, Jjj = 1 if j is odd and 0 
otherwise), (b), Strongly interacting bosons: oscillations are 
completely suppressed and the number of atoms in the mode 
is rather well-defined, although less squeezed than in the Mott 
insulator. {N = L = 6 , U/J = 10, 'y/J = 0.1, Jjj = 1 if j is 
odd and 0 otherwise). (c), Attractive fermionic Hubbard model 
in the strong interaction limit. Measuring only the total pop¬ 
ulation at odd sites = iV-fodd-l-Aj.odd quickly creates singly 
occupied sites, demonstrating measurement-induced break-up 
of fermion pairs {L = 8 , Nj = Aj, = 4:, U/J = 10, 7 /J = 0.1, 
Jjj = 1 if j is odd and 0 otherwise), (d). The same as in (c), 
but with added measurement of the magnetization at odd 
sites Dy = Modd = A^odd — Aj,odd. This protects the dou¬ 
bly occupied sites, thus, demonstrating protection of fermion 
pairs by measurement. The distribution of Aodd vanishes for 
odd numbers, implying that the fermions tunnel only in pairs. 
Simulations for ID lattice. 


l'I'(J.C/.7)) — 


1 + - - - 

A - f47 ^ ^ 

ihj) 


l^n). 


( 10 ) 


The variance of the measurement operator for this state 
is given by 


^AN — 


8J^L 

A2 -b I672 




( 11 ) 


where jz is the filling factor. From the form of the de¬ 
nominator we immediately see that both interaction and 
measurement squeeze with the same quadratic depen¬ 
dence and that the squeezing is always better than in 
the ground state (corresponding to 7 = 0 ) regardless 
of the value of U/J. Also, depending on the ratio of 
7 /A the squeezing can be dominated by measurement 
( 7 /A 1 ) or by interaction ( 7 /A <C 1 ) or both pro- 











cesses can contribute equally (y/C/ » 1). The D = TVodd 
measurement should have a similar dependence on 7 and 
U and be proportional to + 7 ^)”^ since the 7 coeffi¬ 
cient in the perturbative expansion depends on the value 
of (We can see the system transitioning 
into the strong measurement regime in Fig. |^a) as the 
[/-dependence flattens out with increasing measurement 
strength. In typical many-body systems in absence of 
measurement, strong correlations are a product of large 
interactions and the non-interacting limit reduces essen¬ 
tially to a single particle theory. In contrast, here, the 
measurement is another mechanism generating entangle¬ 
ment and strong correlations. This is more evident in 
weak interacting limit, where the measurement-induced 
interaction takes over the standard local interaction (as 
presented in the case of giant oscillations). We will also 
see in the next section how strong measurement leads to 
long-range correlated tunneling events. 

For fermions (Figs. & ,c), the ground state of the at¬ 
tractive Hubbard model in the strong interacting regime 
contains mainly doubly occupied sites (pairs) and empty 
sites. Therefore, in absence of measurement, the fluctua¬ 
tions in the atom population = N^odd + N^^dd 
increase with U/ J while the ones in the magnetization 
by = Modd = N^odd - fV^odd (ci)^) decrease because 
singly occupied sites become more improbable. The mea¬ 
surement induces two different kinds of dynamics using 
the same mode functions since, depending on the light 
polarization, we can address either the total population, 
or both total population and magnetization. In the first 
case (Figs. c), the weak measurement quickly squeezes 
ajj , but destroys the pairs as it does not distinguish be¬ 
tween singly or doubly occupied sites. Figure illus¬ 
trates such a measurement-induced breakup: the initial 
state contains mainly even values of Nodd, corresponding 
to a superposition of empty and doubly occupied lattice 
sites. As time progresses and photons are scattered from 
the atoms, the variance is squeezed and unpaired 
fermions are free to tunnel across the lattice allowing 
odd values for A/^dd- In contrast, probing both density 
and magnetization reduces both their fluctuations and 
increase the lifetime of doubly occupied sites. The result¬ 
ing measurement-induced dynamics is illustrated in Fig¬ 
ure id): Atoms tunnel only in pairs with opposite spin 
as the probability distribution of Aodd contains only even 
values, hence demonstrating the measurement-based pro¬ 
tection of fermion pairs. In both cases, the dynamics of 
the system is not a result of the projective quantum Zeno 
effect (here measurement is weak) but is a manifestation 
of the squeezing of the atomic population and magneti¬ 
zation of the two macroscopically occupied modes. 


V. EMERGENT LONG-RANGE GORRELATED 
TUNNELING 

When 7 ^ J, the photodetections freeze the modes’ 
atom number, and decorrelate the populations of differ¬ 


ent modes. In the quantum Zeno limit of projective 
measurement (7 —> 00 ) tunneling through the bound¬ 
aries between modes would be fully suppressed. However, 
by considering a finite 7 /J we observe additional dy¬ 
namics while the usual atomic tunneling is still strongly 
Zeno-suppressed. In this regime, we observe the follow¬ 
ing effects. First, the evolution between nearest neigh¬ 
bors within each mode is basically unperturbed by the 
measurement process and it is determined by the usual 
tunneling. Importantly, by engineering the light mode 
functions, it is possible to forbid part of this first-order 
dynamics and select which processes participate in the 
quantum Zeno dynamics. Therefore, one can design 
the Zeno subspace of the Hilbert space where the sys¬ 
tem evolves. Second, tunneling between different spatial 
modes is possible only via higher-order, long-range corre¬ 
lated tunneling events that preserve the eigenvalue of b. 
In other words, atoms tunnel across distant sites of the 
lattice via a virtual state and they behave as delocalized 
correlated pairs. 

We can, once more, gain insight into this process from 
the non-Hermitian Hamiltonian. By looking at a second 
order expansion [631164| of the Hamiltonian confined to 
a Zeno subspace of a two-mode measurement, R = 2, for 
U/ J = 0 we obtain the following effective Hamiltonian 


Hz = Po 


ihj) ^ 


blbblh 


Pd 


( 12 ) 


where Pq is the projector into the Zeno subspace, A = 
(Jy, — Jip'b is a constant that depends on the measure¬ 
ment scheme, ip denotes a set of sites belonging to a single 
mode and p is the set’s complement (e.g. odd and even 
sites). 

We can infer a lot of information from this simple ex¬ 
pression. First, we see that first order tunneling will 
only survive between neighboring sites that belong to the 


same mode, because pQb\bjpQ = 0 otherwise. Second, the 
second-order term exists between different modes and oc¬ 
curs at a rate ^ /l- The imaginary prefactor means 

that this tunneling behaves like an exponential decay 
(overdamped oscillations). This picture is consistent with 
the projective limit (7 —>■ 00 ), where the atomic state is 
constrained to an eigenspace of the measurement opera¬ 
tor. 

Crucially, what sets this effect apart from usual many- 
body dynamics with short-range interactions is that first 
order processes are selectively suppressed by the global 
conservation of the measured observable and not by the 
prohibitive energy costs of doubly-occupied sites, as is 
the case in the t — J model [53]. This has profound 
consequences as this is the physical origin of the long- 
range correlated tunneling events represented in (121 by 
the fact that sites j and k can be very distant. This is 
because the projection Pq is not sensitive to individual 
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FIG. 6. Long-range correlated tunneling and entanglement, dynamically induced by strong global measurement in a single 
quantum trajectory. Panels (a),(b) and (c) show different measurement geometries, implying different constraints. Panels 
(1): Schematic representation of the long-range tunneling processes, when standard tunneling between different zones is Zeno- 
suppressed. Panels (2): Evolution of on-site densities; atoms effectively tunnel between disconnected regions due to correlations. 
Panels (3): Entanglement entropy growth between illuminated and non-illuminated regions. Panels (4): Correlations between 
different modes (solid orange line) and within the same mode (dashed green line); atom number Nj {Nni) in illuminated 
(non-illuminated) mode, (a) Atom number in the central region is frozen: The system is divided into three regions and 
correlated tunneling occurs between non-illuminated zones (a.l). Standard dynamics happens within each region, but not 
between them (a.2). Entanglement build up (a.3). Negative correlations between non-illuminated regions (dashed green line) 
and zero correlations between the Ni and Nni modes (solid orange line) (a.4). Initial state: 1 1,1,1,1,1,1,1), y/J = 100, 
Jjj = [0,0,1,1,1,0, 0]. (b) Even sites are illuminated, freezing Neven and Aodd. Long-range tunneling is represented by 

any pair of one blue and one red arrow (b.l). Correlated tunneling occurs between non-neighboring sites without changing 
mode populations (b.2). Entanglement build up (b.3). Negative correlations between edge sites (dashed green line) and zero 
correlations between the modes defined by Neven and Aodd (solid orange line) (b.4). Initial state: |0,1,2,1,0), y/J = 100, 
Jjj = [0,1, 0,1, 0]. (c) Atom number difference between two central sites is frozen. Correlated tunneling leads to exchange of 
long-range atom pairs between illuminated and non-illuminated regions (c.1,2). Entanglement build up (c.3). In contrast to 
previous examples, sites in the same zones (illuminated/ non-illuminated) are positively correlated (dashed green line), while 
atoms in different zones are negatively correlated (solid orange line) (c.4). Initial state: |0, 2, 2, 0), y/J = 100, Jjj = [0, —1,1, 0]. 
ID lattice, U/J — 0, the trajectory is smooth as no photons are detected. 


site occupancies, but instead enforces a fixed value of the 
observable, i.e. a single Zeno subspace. 

Illuminating only the central region of the optical lat¬ 
tice and detecting light in the diffraction maximum, we 
freeze the atom number fViiium imiSI] (Fig. §i). The 


measurement scheme defines two different spatial modes: 
the non-illuminated zones 1 and 3 and the illuminated 
one 2. Figure [^.2 illustrates the evolution of the mean 
density at each lattice site: typical dynamics occurs 
within each region but the standard tunneling between 
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different zones is snppressed. Importantly, processes that 
do not change -/Viiium are still possible since an atom from 
1 can tnnnel to 2, if simultaneonsly one atom tnnnels 
from 2 to 3. Thns, effective long-range tnnneling between 
two spatially disconnected zones 1 and 3 happens dne to 
two-step processes 1—>-2—7>3or3—?'2—J^l. These tran¬ 
sitions are responsible for the negative (anti-)correlations 
{ 6 N 1 SN 3 ) = {N 1 N 3 ) — {Ni){N 3 ) showing that an atom 
disappearing in 1 appears in 3, while there are no nnmber 
correlations between illuminated and non-illuminated re¬ 
gions, {{ 6 N 1 + SN 3 ) 6 N 2 ) = 0 (Fig. 1^.4). In contrast to 
fully-projective measurement, the intermediate (virtual) 
step in the correlated tunneling process builds long-range 
entanglement between illuminated and non-illuminated 
regions (Fig. §r.3). 


To make correlated tunneling visible even in the mean 
atom number, we suppress the standard Bose-Hubbard 
dynamics by illuminating only the even sites of the lat¬ 
tice (Fig. Even if this measurement scheme freezes 
both iVeven and iVoddj atoms can slowly tunnel between 
the odd sites of the lattice, despite them being spatially 
disconnected. This atom exchange spreads correlations 
between non-neighboring lattice sites on a time scale 
~ 7 /J^. The schematic explanation of long-range cor¬ 
related tunneling is presented in Fig. .1: the atoms 
can tunnel only in pairs to assure the globally conserved 
values of A^even and A^oddi such that one correlated tun¬ 
neling event is represented by a pair of one red and one 
blue arrow. Importantly, this scheme is fully applicable 
for a lattice with large atom and site numbers, well be¬ 
yond the numerical example in Fig. .1, because as we 
can see in (121 it is the geometry of quantum measure¬ 
ment that assures this mode structure (in this example, 
two modes at odd and even sites) and therefore under¬ 
lying pairwise global tunneling. Such long-range correla¬ 
tions and long-range entanglement (even in a ID system) 
can develop essentially because of the coherent global ad¬ 
dressing. In contrast, the local uncorrelated probing of 
individual sites would decreases the probability of indi¬ 
vidual tunneling events being correlated. 


This global pair tunneling may play a role of a building 
block for more complicated many-body effects. For ex¬ 
ample, a pair tunneling between the neighbouring sites 
has been recently shown to play important role in the 
formation of new quantum phases, e.g., pair superffuid 
|65| and lead to formulation of extended Bose-Hubbard 
models IM|. The search for novel mechanisms providing 
long-range interactions is crucial in many-body physics. 
One of the standard candidates is the dipole-dipole inter¬ 
action in, e.g., dipolar molecules, where the mentioned 
pair tunneling between even neighboring sites is already 
considered to be long-range |65l 166] . In this context, our 
work suggests a fundamentally different mechanism origi¬ 
nating from quantum optics: the backaction of global and 
spatially structured measurement, which as we prove can 
successfully compete with other short-range processes in 
many-body systems. This opens promising opportunities 
for future research. 


The scheme in Fig. .1 can help to design a nonlocal 
reservoir for the tunneling (or “decay”) of atoms from one 
region to another. For example, if the atoms are placed 
only at odd sites, according to © their tunnelling is 
suppressed since the multi-tunneling event must be suc¬ 
cessive, i.e. an atom tunnelling into a different mode, </?', 
must then also tunnel back into its original mode, ip. If, 
however, one adds some atoms to even sites (even if they 
are far from the initial atoms), the correlated tunneling 
events become allowed and their rate can be tuned by 
the number of added atoms. This resembles the repul¬ 
sively bound pairs created by local interactions ISZIISH]. 
In contrast, here the atom pairs are long-range correlated 
due to the global measurement. Additionally, these long- 
range correlations are a consequence of the dynamics be¬ 
ing constrained to a Zeno subspace: the virtual processes 
allowed by the measurement entangle the spatial modes 
nonlocally. Since the measurement only reveals the total 
number of atoms in the illuminated sites, but not their 
exact distribution, these multi-tunelling events cause the 
build-up of long range entanglement. This is in striking 
contrast to the entanglement caused by local processes 
which can be very confined, especially in ID where it 
is typically short range. This makes numerical calcula¬ 
tions of our system for large atom numbers really diffi¬ 
cult, since well-known methods such as Density Matrix 
Renormalization Group and Matrix Product States (SS] 
(which are successful for short-range interactions) rely on 
the limited extent of entanglement. 


These types of configurations open intriguing oppor¬ 
tunities for quantum engineering of system - bath inter¬ 
actions: here both the system and reservoir, represented 
by different modes, are many-body systems with internal 
long-range entanglement. 


The negative number correlations are typical for sys¬ 
tems with constraints (superselection rules) such as fixed 
atom number. The effective dynamics due to our global, 
but spatially structured, measurement introduces more 
general constraints to the evolution of the system. For 
example, in Fig. we show the generation of positive 
number correlations (shown in Fig. |^.4) by freezing the 
atom number difference between the sites (A^odd ~ A^evem 
by measuring at the diffraction minimum). Thus, atoms 
can only enter or leave this region in pairs, which again 
is possible due to correlated tunneling (Figs. -2) and 
manifests positive correlations. Since this corresponds 
to a no photon trajectory, we can solve exactly for a 
small system and for two atoms, these correlations grow 
as {SfiiSn^) « [1 — sech^(4 J^t/7)]/4. As in the previous 
example, two edge modes in Fig.|^ can be considered as a 
nonlocal reservoir for two central sites, where a constraint 
is applied. Note that, using more modes, the design of 
higher-order multi-tunneling events is possible. 
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VI. SUMMARY 

We proved that the quantum backaction of a global 
measurement can efficiently compete with standard lo¬ 
cal processes in strongly correlated systems. This intro¬ 
duces a physically novel source of competition in research 
on quantum many-body systems. The competition be¬ 
comes efficient due to the ability to spatially structure 
the global measurement at a microscopic scale compa¬ 
rable to the lattice period, without the need for single 
site addressing. The extreme tunability of the setup 
we considered allows us to vary the spatial profile of 
the measurement operator, effectively tailoring the long- 
range entanglement and long-range correlations present 
in the system. The competition between the global back- 
action and usual atomic dynamics leads to the production 
of spatially multimode macroscopic superpositions which 
exhibit large-scale oscillatory dynamics and could be used 
for quantum information and metrology. Such dynami¬ 
cal states show spatial density-density correlations with 
nontrivial periods and long-range coherence, thus hav¬ 
ing supersolid properties, but as an essentially dynami¬ 
cal version. We showed the possibility of measurement- 
induced break-up and protection of strongly interacting 
fermion pairs. In the strong measurement regime, the 
usual nearest-neighbour tunnelling is suppressed but the 
atoms can still tunnel across the lattice because of cor¬ 
related tunnelling. Such globally paired tunneling due 
to a fundamentally novel phenomenon can enrich physics 
of long-range correlated systems beyond relatively short- 
range interactions expected from standard dipole-dipole 
interactions [S5l [55] . These nonlocal high-order processes 
entangle regions of the optical lattice that are discon¬ 
nected by the measurement. Using different detection 
schemes, we showed how to tailor density-density corre¬ 
lations between distant lattice sites. Quantum optical 
engineering of nonlocal coupling to environment, com¬ 
bined with quantum measurement, can allow the design 
of nontrivial system-bath interactions, enabling new links 
to quantum simulations cojin] and thermodynamics m 
and extend these directions to the field of non-Hermitian 
quantum mechanics, where quantum optical setups are 
particularly promising m- Importantly, both systems 
and baths, designed by our method, can be strongly cor¬ 
related systems with internal long-range entanglement. 

Our predictions can be tested using both macroscopic 
measurements [25H271I74] as well as novel methods based 
on single-site resolution [TSHZHj- A pathway to realize 
this is to combine several recent experimental break¬ 
throughs: a BEG was trapped in a cavity, but without 
a lattice p5[ - [27| : detection of light scattered from truly 
ultracold atoms in optical lattices was performed, but 
without a cavity [HI |75], and very recently an optical 
lattice has been obtained in a cavity |28l [55] . Further¬ 
more, the single-atom and multi-particle quantum Zeno 
effect was observed by light scattering 17H1|75|. 

A source of decoherence affecting an experimental real¬ 
isation of our setup is photon miscounts due to imperfect 


detectors. This is more important in the weak measure¬ 
ment regime but its effect strongly depends on the de¬ 
tection scheme (multicomponent states are more fragile 
than single component ones). For strong measurement, 
the exact number of photons detected is not important as 
the state of the system is determined by the photon emis¬ 
sion rate. Atom losses caused by spontaneous emission 
and heating are another source of error that increase the 
uncertainty of the atomic state. Based on off-resonant 
scattering and thus being non-sensitive to a detailed level 
structure, our approach can be applied to other arrays 
of natural or artificial quantum objects: molecules (in¬ 
cluding biological ones) |50|, ions m, atoms in multiple 
cavities [82], semiconductor [83] or superconducting [84] 
qubits. 
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APPENDIX A: ATOMIC HAMILTONIANS AND 
LIGHT-MATTER COUPLING 

We describe atomic dynamics by the (Bose-) Hubbard 
Hamiltonian. For the bosonic case, 

Ho = -hj'^b'^-bi + ^'^ni{ni-l), (13) 

(id) i 

while for the fermionic case 

iJo = -hJ ^ ^ (14) 

<T=T,f (id) i 

where h and are respectively the bosonic and fermionic 
annihilation operators, h is the atom number operator, 
and U and J are the on-site interaction and tunneling 
coefficients. The Hamiltonian of the light-matter system 

is [g 

H = Hq+'^ hLOia\ai + h'^^ima]amFim (15) 

I l,m 

where ai are the photon annihilation operators for the 
light modes with frequencies w;, film = gigm/^ai gi are 
the atom-light coupling constants, and Aq = ujp — uia is 
the probe-atom detuning. The operator Fim = Dim + 
Bim couples atomic operators to the light fields: 

Dim = Bim = Y. JlTblbj, (16) 

* (*J> 

At = uX(r)M^(r)u>(r - r^) dr. (17) 
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Fim originates from the overlaps between the light mode 
functions ui(r) and density operator n{r) = ^^(r)4'(r), 
after the matter-field operator is expressed via Wannier 
functions: 'l'(r) = 'Y^^hiw{T — r^). Dim sums the den¬ 


sity contributions hi, while Bim sums the matter-field 
interference terms. The light-atom coupling via opera¬ 
tors assures the dependence of light on the atomic quan¬ 
tum state. These equations can be extended to fermionic 
atoms introducing an additional index for the polariza¬ 
tion of light modes. 

From (14) we can compute the Heisenberg equation for 
light operators in the stationary limit [5]. Specifically, we 
consider two light modes: a coherent probe beam ag and 
the scattered light oi, which is enhanced by a cavity with 
decay rate k. This allows us to express the cavity mode 
annihilation operator as ai = CFiq, where 


iflioao 
iAp — k’ 


(18) 


Ap = ajg — LUi is the probe-cavity detuning, and C is the 
cavity analog of the Rayleigh scattering coefficient in free 
space |83| . In the main text, we drop the subscript in oi 
and superscripts in J*™. 


ment [30]. It is possible to represent dynamics of a 
system with two different processes: non-Hermitian dy¬ 
namics and quantum jumps. The first one describes 
evolution of the system between two consecutive mea¬ 
surement events, while the second one models the pho- 
todetections. These can be simulated introducing the 
jump operator c = and effective Hamiltonian 

HeS = Hg — ihc'cf2. Starting from the time tg, a ran¬ 
dom number r between 0 and 1 is generated with uniform 
probability, and we solve the effective Schrodinger equa¬ 
tion 


±\^(t)) = (19) 

until the time moment tj such that 
At this time moment, a photon is detected and the jump 
operator is applied to the state of the system, which is 
subsequently normalized: 




c|V'(^j)) 

V'(V'(ij)|clc|V'(ti)) 


( 20 ) 


APPENDIX B: MODELING THE 
MEASUREMENT PROCESS 

Quantum trajectories describe the evolution of a quan¬ 
tum system conditioned on the result of a measure- 


Finally, a new random number is generated and the pro¬ 
cedure is iterated setting tj as a new starting time. This 
technique allows us to simulate a single quantum trajec¬ 
tory, thus modeling the result of a single experimental 
run. 
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